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Abstract 

We present an improved form of the integration technique known as NDIM (Negative Dimensional 
Integration Method), which is a powerful tool in the analytical evaluation of Feynman diagrams. Using 
this technique we study & (f>'^ (B 4>* theory in D = 4 — 2e dimensions, considering generic topologies 
of L loops and E independent external momenta, and where the propagator powers are arbitrary. 
The method transforms the Schwinger parametric integral associated to the diagram into a multiple 
series expansion, whose main characteristic is that the argument contains several Kronecker deltas 
which appear naturally in the application of the method, and which we call diagram presolution. The 
optimization we present here consists in a procedure that minimizes the series multiplicity, through 
appropriate factorizations in the multinomials that appear in the parametric integral, and which 
maximizes the number of Kronecker deltas that are generated in the process. The solutions are 
presented in terms of generalized hypergeometric functions, obtained once the Kronecker deltas have 
been used in the series. Although the technique is general, we apply it to cases in which there are 
2 or 3 different energy scales (masses or kinematic variables associated to the external momenta), 
obtaining solutions in terms of a finite sum of generalized hypergeometric series de 1 and 2 variables 
respectively, each of them expressible as ratios between the different energy scales that characterize the 
topology. The main result is a method capable of solving Feynman integrals, expressing the solutions 
as hypergeometric series of multiplicity {n — 1), where n is the number of energy scales present in the 
diagram. 
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1 Introduction 



In Quantum Field Theory the permanent contrast between experimental measurements and theo- 
retical models has been possible due to the development of novel and powerful analytical and numerical 
techniques in perturbative calculations. The fundamental problem that arises in perturbation theory is 
the actual calculation of the loop integrals associated to the Feynman diagrams, whose solution is spe- 
cially difficult since these integrals contain in general both ultraviolet (UV) and infrared (IR) divergences. 
Using the dimensional regularization scheme, which extends the dimensionality of space-time by adding a 
fractional piece (D = 4 — 2e), it is possible to know the behavior of such divergences in terms of Laurent 
expansions with respect to the dimensional regulator e when it tends to zero. On the other hand, the 
structure of the integral associated to the diagram gets increasingly more complicated when the number of 
external lines, loops or energy scales is increased, and therefore finding an analytical solution is extremely 
difficult. Of the many different techniques [1] that have been developed in order to evaluate diagrams, 
we can mention: integration by parts (IBP), contour parametric integrations in terms of a Mellin-Barnes 
representation of the diagram, the differential equations method (DEM), etc. 

In the context of the present article, and as comparison, particular mention has to be made to the 
Mellin-Barnes integral representation of the diagram, which is a very powerful technique for finding solu- 
tions in cases of arbitrary propagator powers. These powers can typically occur when tensorial structures 
are present, when in the reduction process generate scalar integrals with different powers of propaga- 
tors, or when there are subtopologies associated to massless virtual particles contained in a propagator, 
which adds to the corrected propagator a fractional piece proportional to the dimensional regulator e. 
Something similar happens when the Integration by Parts (IBP) technique is used. The Mellin-Barnes 
integral representation consists in transforming the momentum space representation of the diagram into 
a parameter space representation, written in terms of multiple contour integrals. As solution of these 
contour integrals one obtains generalized hypergeometric series, whose arguments can be T' or ratios of 
the different energy scales that are present in the diagram, and which give information about the different 
kinematical regions, specified in a natural way by the convergence conditions of these functions. The most 
interesting aspect of these type of solutions is that they can be analyzed and expanded in terms of the 
dimensional regulator e in the kinematical region of interest. For this purpose there are known techniques, 
algorithms and calculation packages [5] , [3] , [J] . 

The specific technique that is used here was suggested originally by Halliday and Ricotta [Sj , and it is 
known as NDIM (Negative Dimensional Integration Method) since it performs an analytical continuation 
in the dimension D to negative values. This can be done due to the fact that the loop integrals have the 
property of being analytical functions in arbitrary dimension. In actual practice this technique represents 
the diagram in terms of a multiple series whose argument contains a certain number of Kronecker deltas, 
and the solutions emerge naturally from evaluating the sums in the different forms that these Kronecker 
deltas allow. The final result can be always expressed in terms of generalized hypergeometric series. This 
method has been developed and used by several authors, among them A. Suzuki, A. Schmidt, E. Santos, 
C. Anastasiou, E. Glover and C. Oleari, in many applications to essentially one-loop diagrams. 

Since the results obtained in both the Mellin-Barnes representation of diagrams and in this technique 
are generalized hypergeometric functions, it is possible to compare them. In these series solutions the 
arguments are ratios of the different energy scales that appear in the diagram, and the parameters that 
characterize them are linear combinations of the propagator powers and the dimension D. In such appli- 
cations it was shown that at one-loop level [H [71 [H [H [10] this technique is comparable in terms of the 
complexity of the solutions, with the one coming from the Mellin-Barnes representation. Nevertheless, in 
two or more loop cases [TH [T3I [HI [ISl [HI [H] , the conclusion is that the generalization of the formalism 
is not adequate since the solutions get to be complicated in two essential aspects. First, the number of 
terms or contributions of the obtained hypergeometric functions is very high, being for some two-loop 
diagrams typically of the order of thousands; and second, the complexity of each of them increases, due to 
the fact that the multiplicity of summations of each hypergeometric series is also quite high and goes over 
the number needed to represent the solutions with respect to the number of kinematical variables that the 
diagram possesses. Both of these aspects make that any analytical study of the solutions to the L-loop 
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case be very complex, and therefore it becomes impossible to identify each of the contributions associated 
to a specific kinematical region. Nevertheless, this technique has an important advantage, since it is much 
simpler to apply that the Mellin-Barnes representation, because it changes the process of solving multiple 
contour integrals into solving a linear system of equations, a fact which is very important when one is 
trying to find analytical solutions. 

We have verified that in the form presented by the above mentioned authors, the application of this 
method to the general L-loop case is indeed very cumbersome. This can be observed in the L — linear 
and (L + 1) — linear multinomials that are present in the Schwinger parametric integral, which contains 
a number of terms that increases rapidly when the number of loops or independent external momenta 
increases. Therefore the number of multiple expansions grows with respect to the number of constraints 
or Kronecker deltas which can be obtained from the mathematical structure of the diagram and which 
are an essential part of the method. In the present work we have deduced the method in a different way, 
and in the process found analytical solutions to L — loop diagrams which moreover can contain massive 
propagators, extending significantly in this manner the number the diagrams that can be treated with 
this technique. In our formalism it is not necessary to work with the momentum representation of the 
graph. We start with its Schwinger parametric representation, using for this purpose the mathematical 
structure of a generic L loops, N propagators and E independent external lines diagram, shown in Ref. 
[IH], and which provides a quick and direct way of obtaining the parametric representation integrand 
without explicitly solving the momenta integrals. We concluded that this technique is quite powerful and 
has significant advantages in the evaluation of an infinite set of diagrams belonging to certain topological 
classes, which is something that will be shown in the specific examples that are going to be presented in 
this work. 

Our emphasis in this work will be mainly in finding appropriate procedures for obtaining the general 
solutions associated to a diagram, instead of analyzing in more detail each one of them, although we 
do it in cases in which they are naturally related to the convergence relations of the hypergeometric 
series, which in turn define the kinematical regions where these series are valid. The presentation of this 
work is as follows: section (//) contains a description of the algebraic components of NDIM that are 
required for Feynman diagram evaluation, with emphasis on the basic formulae of the method, which are 
related to delta Kronecker generation, topology related series expansions, and the reasons for naming this 
series multiregion expansion. We also explain here how to obtain in a simple way the explicit Schwinger 
parametric representation in the different energy scales of a general topology [18j . Finally we describe 
an algorithm that sums up the above in a systematic sequence of steps for solving Feynman 's integrals 
and thus obtaining the multiregion expansion of the diagram, or which could be called presolution of the 
graph. In section (///) we have considered three examples where the technique is applied systematically, 
starting from a simple known case until a complex one with until now unknown solutions. The first is 
the Bubble diagram [E ^ \,N — 2,L = 1), with equal mass (to) propagators. This is a basic example, 
whose solutions are well known [191 EHl EI] • The following example corresponds to the diagram called 
CBox [E = i^N — b,L = 2), which although has known solutions [11], in the on-shell massless case this 
is the first time in which they have been obtained using this method. Increasing the complexity, we solve 
this diagram for the case of three massive propagators, which by itself is a new result. Finally we find 
the solution for a four-loop propagator case {E = 1, N — 8, L — i), where the idea is to show that this 
technique is very efficient and effective when one factorizes and expands systematically the multinomials 
of the parametric representation. Initially the massless propagator case is solved, and then with two 
massive propagators. The massless case is quite simple, since it is a diagram which is reducible loop by 
loop; but the second case is much more difficult, and the solution is in terms of hypergeometric series of 
the type gFg, in the variables {p^/Arn^) or its reciprocal [Am'^/p'^), according to the kinematical region 
of interest. This result is relevant given the difficulty of obtaining it using any other method in the case 
of arbitrary propagator powers. In section {IV) we discuss some conceptual aspects of the integration 
method and about the complexity of the solutions when the number of loops L increases. We have added 
a section {V) as an extension of (III), in which we present several topologies that contain two or three 
different energy scales, and where we only show the presolution in each case. Finally in the appendices 
we include a summary of definitions and properties of one and two variable hypergeometric functions. 
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which appear frequently in the solutions of Feynman integrals; and also a more detailed description of the 
equations that support the integration method here presented, with some of its properties. 

The most important aspect of this work is that we have been able to associate this technique with a 
certain family of L-loops topologies, with respect to which it is applied very easily and directly, with the 
possibility of solving a large number of Feynman diagrams, which can even contain massive propagators 
with arbitrary powers. This allows to extend considerably the known solutions of Feynman diagrams. 



2 Mathematical foundations 

In order to understand the method, we describe now its main components: 

- Obtaining Schwinger's parametric representation.. 

- Foundations of the integration method NDIM. 

- General algorithm. 

2.1 The parametric representation 

Let us consider a generic topology G that represents a Feynman diagram in a scalar theory, and 
suppose that this graph is composed of N propagators or internal lines, L loops (associated to independent 
internal momenta q — {qi, ...,qi,}^ and E independent external momenta p = {pi, ...,pe}- Each propagator 
or internal line is characterized by an arbitrary and in general different mass, to = {rrii, to^v}- 

Using the prescription of dimensional regularization we can write the momentum integral expression 
that represents the diagram in Z? = 4 — 2e dimensional Minkowski space as: 

r-r( \^ rf'^gL 1 1 f.. 

In this expression the symbol Bj represents the momentum of the j propagator or internal line, which in 
general depends on a linear combination of external {p} and internal {g} momenta: Bj — Bj{q,p). We 
also define i/ = {vi, i^n} as the set of powers of the propagators, which in general can take arbitrary 
values. After introducing Schwinger's representation, it is possible to solve the momenta integrals in terms 
of Gaussian integrals. The result is Schwinger's parametric representation of the diagram associated to 
equation ([T]) , which in the general case is given by the expression: 



(-1)-^ r ^^pfe=i^^"^0^^p("c7 

G=^N^ dl^ ^ ^ (2) 

We have introduced for simplicity the following notation d~x — dxi...dxN Il^i where U and F 

are L — linear and (L + 1) — linear multinomials respectively, defined in terms of determinants [181: 

E 



^ . (L+1) 

(3) 



t^iL.,pf+^j:t^iL.,p^-p. 

i—l i—1 j>i 



where A^*^^^' is the determinant defined by the equation: 
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Ak+1) 



Ma 



Mik M: 



Mik Mij 



(4) 



The matrix symmetric M has dimension {L + E) x {L + E), and is called Initial Parameters Matrix. It 
is possible to easily built it when ([1]) is parametrized, and the internal products of loop momenta and 
external momenta are expanded, with coefficients which correspond to the elements of the matrix My. 
For a better understanding of this process let us define the momentum: 



if L>j>l, 



Pj-L if {L + E)> J > L, 



(5) 



with the (L + E)-vector Q = [Qi Q2 ... Q(^l+e)Y ■ Using this definition the following matrix structure is 
generated, as a previous step to the loop momenta integration: 



L+E L+E 

i=i j=i 



(6) 



with which we finally identify the symmetric matrix M. Let us develop briefly a simple example, for the 
massless one loop diagram shown in Fig. 1. 




Figure 1: Bubble diagram. 
The corresponding 1-loop integral is given by the expression: 

f d^q 1 

J inD/^q^q+pf 

We then apply Schwinger's parametrization, and then obtain the following equation: 



G = J .^0% J dxi exp (-xi g^) J dx2 exp (^-X2 {q + pf^ 



(7) 



(8) 





or equivalently, expanding the squares in the exponents and factorizing in terms of internal and external 
momenta, it follows that: 



G = dxidx2 



d^q 



exp [- [(a;i + X2) q^ + 2x2 Q-P + X2 p'^]] 



(9) 



Now the matrix of parameters is obtained directly from the exponent by simple observation, that is: 

M.(-;--). ,10) 

where it has been used that Qi = q and Q2 = p. Later on other examples are going to be presented. For 
more details see Ref. [T51. 
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2.2 Foundations of the integration method NDIM 

2.2.1 Notation and essential formulae 

The exponential structure that the integral of parameters ((21) presents allows that the integration 
method NDIM be developed and sustained by two formulae that can be deduced starting from the following 
integral expression: 



C30 

1 1 /" 



(11) 



where the quantities A and (3 are in general complex. Using this equality one can justify the essential 
point of this method, and which is related to the equivalence (from an operational point of view) of the 
integration sign with a quantity which is proportional to a Kronecker delta. In order to show this we 
expand the integrand of equation 

oo 

n p 

Now the corresponding evaluation of the integral will not be done in the usual way, but we define an 
identity which fulfills the equality. This happens after we define the following operational relation: 

oo 

j dx xP+--^ ^ r (/5) h+nfi- (13) 



For convenience we introduce the following notation in order to write the operational equivalent of the 
integral sign: 

= + (14) 

where the parenthesis (■) has implicit the constraint associated to the Kronecker delta and which further- 
more satisfies the following property (see appendix): 

(« + /?) EE r(~c.)^|^±^<5„+^,o, (15) 



or in a simplified manner: 



(a + /3) = Wo, (16) 



where u is an arbitrary index or parameter and where we have also defined the factor: 

(-1)" 



(17) 



On the other hand, by applying successively equation (jll[) and then equation ^\A\ to an arbitrary multi- 
nomial, we can deduce the second fundamental formula on which the method is based. This expresses 
the fact that a multinomial of a terms can be written as a multiregion expansion, in such a way that it 
simultaneously contains all limiting possible cases with respect to the relation between the different terms 
present in the multinomial. The series so described can be written as: 

(^1 + - + A.f' <^n....^. AT...A-'^^^^^^^±-±^, (18) 
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where the terms Ai {i — 1, ...,a) and the exponent v are quantities that can take arbitrary values. The 
expression contains the a expansions that can be obtained and which have the general form: 



' ni=0 n„=0 \ * / \ 

where the definition of the factor in (fT7|) has been generalized: 




(19) 



In this case each of the expansions (fTO)) that can be obtained starting from (fT5|) correspond to multivariable 
generalized hypergeometric series of multiplicity pt— (a — \\ and each one of them contains the limiting 

cases or region where it is fulfilled that ^^<l^Vj^i. 

2.2.2 General form of the multiregion expansion (MRE) of a diagram and its solution 

Once the parametric representation of a diagram ^ has been obtained, the following step is to 
evaluate it, and for this purpose it is necessary to make an expansion of the integrand starting from the 
exponential functions if they are present, and then continuing with a multiregion expansion of all the 
multinomials that the process generates according to formula (llSp . The expansion process ends when one 
finally obtains only one term, which is a product of all the Schwinger parameters. Now it only remains 
to replace all the integrals according to the formula (fT4|) or its equivalent (•). We have obtained in this 
manner the presolution or multiregion expansion of the Feynman integral considered in ([1]) . In the case in 
which we consider a generic topology G, characterized by M different mass scales, iV propagators and P 
Lorentz invariants associated to internal products of the external independent momenta, then the general 
form of the multiregion expansion is given by the following expression: 

G = i-i)-'^ E ^rn..n. uiQ'.r n (-™ r n %7^^ n ^1^, (21) 

where it is possible to identify the following quantities: 

(T Corresponds to the multiplicity or number of summations that conform the multiregion expan- 
sion, or in this case the presolution of the diagram G. 

=^ Corresponds to a kincmatical Lorentz invariant, which is a quadratic form of the independent 
external momenta. 

aj, Pj, "fj =4> Correspond to linear combinations of the indexes {ni, .., n^}, with the exception of /S^, 
which includes a dependence on the dimension D: 

/3i = y + ni + ... + 7ip. (22) 

The coefficients of the sum indexes {rii} in the linear combinations aj and 7^ are (+1) and in the case of 
f3j the indexes have coefficients (-1), except for P^. 

N =>Numbcr of propagators or equivalently number of parametric integrations that the method 
transforms into N Kronccker deltas. 

K =>Number related to the total number of multiregion expansions performed over the integrand of 
the parametric representation, which in turn generates K constraints or equivalently K Kronccker deltas. 
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A particular case corresponds to vacuum fluctuation diagrams, in which case the presolution (|2ip is 
written in the following form: 



G = i-i)-'^ v 0„„..,„. n(-™?r- n ^^^^^ n ^^i^, (23) 



where = y- 

In order to find the solutions, it is necessary to evaluate the sums using for this purpose the existing 
constraints among the indexes of the sum, represented by the 6 ~ {N + K) Kronecker deltas; but as it 
can be seen there are several ways to do this evaluation. In fact the number of different forms to evaluate 
the presolution of G using the Kronecker deltas is given by the combinatorial formula: 

Each one of these forms of summing using the Kronecker deltas in G, generates in turn as a result a 
multiple series, which corresponds to a generalized hypergeometric function, whose multiplicity is given 
by: 

/i = (cr - <5) . (25) 

In general it is not always possible to use the 6 deltas in order to evaluate the corresponding S sums, 
since this will depend on the combination of indexes with respect to which the sum should be done. If 
this happens, these cases are simply not considered as contributions to the solution and therefore are 
discarded. In those cases in which we do have a relevant expansion contribution, this corresponds to an 
serie representation of the set of kinematical variables that are present in the problem, in the form of 
ratios of the different energy scales that appear in the topology. Since the multiregion expansion (|2ip 
contains all the limits simultaneously, all the solutions that have been found are related between them by 
analytical continuation, which is realized implicitly by the integration method NDIM. 

In practical terms the idea of the method is to generate finally an expansion that represents the 
diagram G, the Multiregion Expansion, characterized by multiplicity cr, in combination with a number S 
of Kronecker deltas. Using this it is possible to find the solutions in terms of generalized hypergeometric 
functions of multiplicity fi. In this work we study in detail only the topologies where the fi variables of 
the different solutions correspond to ratios of the different {P + M) energy scales present in the topology, 
that is when the following condition if fulfilled: 

H={P + M-1). (26) 

2.2.3 Minimizing expansions and maiximizing the number of constraints (Kronecker deltas) 

This is an essential point in the improvement of the integration technique NDIM, since it explains 
the procedure used to find the optimal multiregion expansion of the diagram G. The idea is quite simple 
and it is related to the form in which the expansion of the multinomials F and U of the parametric 
representation ^ is done. The technique consists in the factorization of all the multinomials which are 
repeated two or more times in the integrand, without expanding them until they factorizc in one term as 
product of the previous expansions. In the same way the process is repeated for all the submultinomials 
that the process generates and ends when finally one obtains one single monomial which is a product of 
the N Schwinger parameters. Finally applying equation (fT4|) . the integrals get associated to N Kronecker 
deltas. 

The main result that is achieved by such a procedure is the minimization of the number of expansions, 
increasing the number of submultinomials that have to be expanded and since each multiregion expansion 
done over them generates one Kronecker delta, using formula (llSp it is clear that the number of deltas is 
also maximized. 
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Let us see the following examples where it is verified that this systematic expansion indeed generates 
generalized hypergeometric functions of minimized multiplicity /i, and that in general less terms than in 
previous approaches are part of the solution of the diagram. 



Example I Let us consider the following function: 

g = {aia2 + aia^ + a2a4 + asQi)^ , (27) 

and let us find its corresponding multiregion expansion, with and without the factorization previously 
indicated. 



Expansion without factorization : Applying directly formula ^8|) to function g, one easily obtains 
the multiregion expansion of the function: 



9= E -^.n.-.n^ <^^"^«r 



"3 "4 



r(-/3) 



Expansion with factorization : We now reorder the function g, factorizing the multinomial in the 
following way: 

g^laib + a^bf, (28) 

where b = (02 + a^). We have applied the previously stated idea with respect to repeated submultinomials. 
We now expand the binomial in (pS)) . obtaining the following: 

/."i+"2,«i n2 {-P + ni+n2) 
9 — ?*ni,«2 ° "^1 "^4 T(-B) ' ^ ' 

ni,"2 

in an analogous manner the binomial b ~ (a2 + 0.3) is now expanded, which gives as result the multiregion 
expansion of g: 

n- rh .-,rxi„n3„n4„»2 ('/^ + »1 + ^2) ("^i - 7^2 + ^3 + ^4) . . 

■9 - ^^2^^^^n,.-,n, «1 ^2 «3 «4 r(-ni - n2) ' 

The advantage of expanding after using the factorization of repeated submultinomials can be clearly 
appreciated in Table L 



without factorization with factorization 



Multiplicity multiregion series (cr) 4 4 

Kronecker deltas associated to the expansion (6) 1 2 

Multiphcity of resulting series = a — S) 3 2 

Possible resulting expansions (Cg) 4 6 

Relevant resulting expansions 4 4 

(Table I) 

Notice that the multiplicity of the produced hypergeometric series has decreased, although the number 
of relevant contributions, those that correspond to a limiting case of g, has remained the same. 

Example II Let us consider a second example, a function g which is a product of two monomials: 

5 = (ai + 02)" (ai + 02 + as)'^, (31) 
and find its respective multiregion expansion with and without factorization. 
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Expansion without factorization: We expand each multinomial separately, and then after re- 
ordering obtain the following series for g: 



Expansion with factorization : Let us now see what happens if the repeated submultinomials in 
((3T|) are factorized as follows: 

g = + 03]'', (33) 
where h = (ai + 02) . Then expanding the binomial, one obtains the series: 

and then expanding the factor b — (ai +02), which finally gives us the multiregion expansion for g: 

ff= 2. «2 «3 — rF^^) fF^^ ■ (35) 

ni,..,n4 ^ ' \ ' / 

The important fact is that once again the complexity of the resulting series representations gets reduced, 
both in the multiplicity /i and in the number of relevant contributions, as can be seen in Table II. 



without factorization with factorization 



Multiplicity of multiregion series (a) 


5 


4 


Kronecker deltas associated to the expansion (<5) 


2 


2 


Multiphcity of resulting series {fi = a — S) 


3 


2 


Possible resulting expansions (Cf ) 


10 


6 


Relevant resulting expansions 


6 


5 



(Table II) 

Although we have presented very simple examples, they illustrate the advantages of first doing the factor- 
ization and then expanding systematically a multinomial. In these examples we verify that the multiplicity 
of each hypergeometric function finally obtained decreases when using the factorization procedure of re- 
peated submultinomials, and also the number of relevant expansions obtained from the multiregion series 
is reduced as well. 

For the particular case of Feynman integrals, the equivalent Schwinger parametric representation is 
composed of two multinomials, F and U, and the factorization and expansion process shown in the 
examples is directly applicable to them. 

2.3 The algorithm 

Here we present an algorithm for finding the solution of an arbitrary Feynman diagram: 

1. Find Schwinger 's parametric representation ([2|) of the Feynman diagram G, characterized by L loops, 
N propagators, E independent external momenta and M different masses. In general {M < N). 

2. Minimize the number of terms of the multinomial reordering such that: 

^ = E E ^iL+n p^-pj = jl /.(^) Ql (36) 

i=l j = l 3 = 1 
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where the factors fj(~x) correspond to multinomials which depend only on Schwinger's parameters, 
and where is a kinematical invariant which depends on the independent external momenta. Thus 
the integrand of the parametric representation acquires the following structure: 

t 2\ ( 2 ^ ( fiQi\ ( fpQ% 
exp (yimf) ... exp [yMm%.j) exp — ... exp ' 



U J "V U 



urn • 

Being {yi, ?/m} C {xi, Xn}. In principle, assuming that all masses are different and non- 
vanishing, (M + P) independent (between them) expansions are obtained, only associated to these 
exponential functions. 

3. The following step is finding repeated multinomials susceptible to factorize, both in the functions 
fj{~x ) as well as in the multinomial U. For the equal masses case, before expanding the exponential 
which contains them, it is convenient to factorize the terms that build a multinomial already existent 
in the previously factorized integrand, and then one proceeds to expand such exponential. 

4. Expand the multinomials until finally a single product of Schwinger's parameters is obtained. Each 
performed expansion will be associated to a Kronecker delta. 

5. After doing all the expansions, finally it remains to replace the integral signs by its equivalent (■). 
By doing that we add N additional Kronecker deltas. The result is the presolution of the diagram 
G or its equivalent multiregion expansion (|2ip. 

6. In order to find the contribution or serie representation associated to a particular combination of 
the /i free indexes, the most appropriate procedure is to solve the linear system corresponding to the 
constraints among the indexes, obtained from the Kronecker deltas, and assume that such indexes 
in these equations correspond to independent free variables. Thus we obtain a set of solutions for 
the indexes that are not free, in terms of the fi free indexes, of the parameters {vi, vn} and of the 
dimension D. Each of these combinations constitute a generalized hypergeometric function, whose 
series representation has multiplicity /i. Not all the combinations of free indexes generate a solution, 
in which case the associated linear system simply has no solution or equivalently the Kronecker 
deltas cannot eliminate the 6 sums in the non-free indexes, so the quantity Cg is an upper bound 
with respect to the number of possible hypergeometric contributions that are present in the solution 
of the diagram G. 

7. Repeat the previous process for the Cg forms of combining the fi indexes. Thus we obtain at 
most Gg hypergeometric series, which are classified according to kinematical region of interest. The 
solution in each kinematical region corresponds to the algebraic sum of all the contributions that 
have the same kinematical argument or variable. The final result is the evaluation of the diagram 
G in terms of all its series solutions. 



3 Applications 

In order to show explicitly the integration technique, let us consider three applications. The first 
corresponds to the evaluation of the one-loop massive propagator. Through this simple and known problem 
we present the formalism and proposed notation. In the second example, the CBox diagram, it is already 
possible to notice the efficiency of an adequate factorization. In this case the solution is also known 
and has been already obtained using a different method to the one proposed here, which is useful in 
order to compare the simplicity of the NDIM with respect to other methods. We also generalize this 
problem to cases in which the solution has not been found until now. And the last example is a four-loop 
diagram which has associated a very complex integral given the number of terms that are present in the 
polynomials F and U, and which nevertheless with an optimal factorization becomes possible to be solved. 
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even considering one mass scale in the graph. Is here where the integration technique demonstrates aU its 
power, obtaining a completely new result and showing the simplicity of the method. 



3.1 Example I : Massive Bubble diagram 

This is a case in which the solutions are well known [191 HOI E] , and which we will use in order to 
introduce the integration method: 




Figure 2: Labelled Bubble diagram. 
The integral representation of this diagram (Fig. 2) in momentum space is given by the equation: 



d% 1 

Let us consider the case of equal mass propagators, mi = m2 = m, and then according to equation ^ 
the corresponding Schwinger parametric representation is given by the expression: 

D °° exp (- (xi +2:2) (-TO^)) exp ( ^iifH_p2 

G^±^^ dH! , ^ ^^+^^ - , (39) 



where we readily identify the multinomials of the representation, that is ?7 {xi + X2) and F — X1X2 Pi- 
Then the only repeated multinomial coincides with U. For finding the multiregion expansion we expand 
first the exponentials, therefore obtaining the series: 

(-1) 

1 

Now U is expanded: 



G = \ y <\>^^ „^ {/T [-^^T / — • (40) 



1 ^ a;"^X2* 

(xi + X2) " „3,„4 i 2 + "-1 "-2j 

which replaced in (|40p and doing a separation of integration variables gives us: 

G = E ^?¥^^^^Ai / dx, .-+-+--1 I (42) 

Y{,=i^i^j) nt^^n, r(^+ni-n2) J J 

Then using equation (fT4l) the integrals are transformed in its equivalent (•), which finally allows us to 
obtain the multiregion expansion of the diagram G: 

G - E ivir {'-^T ri^h^\y (43) 

Y[j = l^{Vj) nu..,n, r(^+ni-n2) 



12 



where we have defined a notation for constraints {A^}: 

Ai = ( Y + ni - ?i2 + n3 + 714) , 

A2 = {1^1 +ni+ ng) , (44) 
A3 = (1^2 +ni+ Tii) . 

The number of possible representations is given by the combinatorics C| = 4, and the obtained hyper- 
geometric series will have multiplicity one. The constraints in ()44p provide us with the following linear 
system: 

= -y + rii — 712 + na + 714, 

Q = vi+ni + 713, (45) 



= i^2 + '^i + n^. 

Explicitly the combinatorics expresses that we have three equations and four variables (indexes), and 
therefore it is necessary to leave one free or independent index, which can be done in four different ways. 
Let us see this case by case: 



3.1.1 Serie representation when Tti is taken as independent index: 

Before doing any calculation let us define a notation for relating the index that is taken as free and 
the respective hypergeometric representation that it generates. For this purpose let us identify Gj as the 
contribution obtained once we leave index Uj free in the multiregion representation of diagram G. 

Now starting form the multiregion expression (|43p . and using formula (|16p . we can replace the paren- 
thesis (•) conveniently, thus obtaining the series with index of sum ni. 

G4 = E ^ {r>ir i-^T ^^'^f^""-^^^T-\ (46) 

n'=ir(^,)tr r(f+77l-772) ' 

where the dependent indexes (solutions of the system in P5)) ) take the following values: 

n2 = Y - '^1 - ^2 - ni, 

113, = -vi~ni, (47) 



Replacing now in ([^^ . we get: 



D —777 



2W 



-f 1+1^2 



n=ir(^,) 



771 



T{vi + V2 + 2ni) 



(48) 



and we can then reduce the above expression in terms of a hypergeometric function 

D 



2\^--''i+^2 p / l'l+V2-^, Vl-, V2 
2 2 ' 2 



A.) 

Am? J 



where: 



Xi = (-1) 



_z,r(i/i + j/2-f) 
2 — 

T(vi + V2) 



(49) 



(50) 
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3.1.2 Serie representation when 712 is independent index: 

Analogously to the previous case, we now make the index 712 independent, and from an adequate 
writing of the factors {A;}, we get the following expression associated to a free 712: 



Go 



n=ir(-. 



] 1 n2 



r(f +ni-n2) 



This time the solutions for the dependent indexes are given by: 



(51) 



ni = 



n2, 



m = V2 - ^- + n2, 



(52) 



D 



"2, 



which allows to obtain the serie representation when 712 is free in the presolution 



G2=X2 (P?)"""^" 3F2 

where the prefactor X2 given by the identity: 

X2 = (-1) 



V1 + V2 - 2 ^ 2 T' ' 2 



2 



r(^i + ^2-f)r(f -t^i)r(f -7.2) 



r(z/i)r(i.2)r(i? - - i.2) 



Am 



(53) 



3.1.3 Solution for and 774 independent: 

Similar procedures to the previous ones allow us to quickly get the terms G3 and G4 respectively: 



2\ -2 -"2 



3F2 



^1, i + 1 
1 + 1^1 - 1^2, 1 



2 
1^2 



and: 



U2. i + ^ 



1 1 

' 2 



where the following factors have been defined: 



X3-(-l) 

Xa = (-1) 



I + V2 - Vl, 1 + -J - I'l 



r(^2) ' 

-iz r(^i-f) 
r(^i) ■ 



4m^ 



(54) 
(55) 

(56) 
(57) 



3.1.4 Solutions in the different kinematical regions, expressed as sums of terms Gj 



We can distribute the previously found solutions in the two regions where it is possible to expand the 

4to^ 

< 1, where the solution of G is given by the expression: 



kinematical variables, the first located in 



G 



9 

Pi 
Am? 

Pi 



G2 + G3 + G4 



(58) 



and the solution in the region where 



4m^ 



< 1: 
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(59) 



In this way we have evaluated G in terms of hypergeometric functions, which correspond naturally to 
serie representations with respect to the two energy scales present in diagram G. 

3.2 Example II : CBox diagram (On-Shell case) 

The next diagram has known solutions [11], found by a different method. An attempt using NDIM 
can be found in Ref. [55]. Nevertheless, the result presented here is an improvement in the sense that it 
is much simpler and easy to obtain. We do not need to use any simplification tools in order to get the 
final result. 





Figure 3: Labelled CBox diagram. 
For this diagram (Fig. 3) the integral representation in momentum space is given by: 

rf^gi d^q^ 1 



G = 



1 



1 



1 



m 2 ITT 

where the quantities Bi correspond to: 



^ ■ * (Bi)-i (Ba)'^^ (B^)-^ (BiY^ [B^Y^^ 



(60) 



(61) 



Bi = qj - mj + iO, 
B2 = {qi + pi)'^ - ml + iO, 
B3 = {qi +q2+Pi+P2f -ml + iO, 
Bi = (qi +q2+Pi+P2+ p-if -m\ + iO, 
B5 = qj- mi + lO. 

According to the method developed in Ref. [18j . we can find the initial parameters matrix associated to 
the topology, which is: 

/ Xi + X2 + X3 + X4 X3 + X4 X2 + X3+ X4 X3 + X4 X4 \ 

X3 + X4 X3+ X4+ X5 X3 + X4 X3+ X4 X4 

M — X2 + X3 + X4 X3 + X4 X2 + X3+ X4 X3 + X4 X4 ■ (62) 

^3 + X4 X3 + X4 X3 + X4 X3 + X4 X4 

y X4 X4 X4 X4 X4 J 

Starting from this matrix we can find the algebraic components of the parametric representation, evalu- 
ating for this purpose the polynomials U and F. The multinomial (2 — lineal) U will be evaluated using 
the determinant: 
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u 



Xi + X2 + Xs + Xi Xs + Xi 

X-i +X4 X3+X4 + X5 



= X5X1 + X5X2 + X5X3 + X5X4 + XiXs + X1X4 + X2X3 + X2X4, 

and the multinomial (3 — lineal) F will be written as: 

F = 2C1.2 PI-P2 + 2Ci,3 Vl-P?, + 2C2,3 P2-P3, 



(63) 



(64) 



where the coefficients Ci^j are evaluated in terms of subdeterminants of the matrix of parameters, and 
then: 



Ci,2 = 
Ci,3 = 
(^2,3 = 

We then have for F: 



Xi + X2 + + X4 
Xi + X4 
Xi + X4 

X\+X2 + Xz + X4 
Xz + X4 

X4 

Xi + X2 + X3 + X4 
X3 + X4 
X4 



X3 + X4 
X3 + X4 + X5 
X3 + X4 

X?, + X4 
X3+X4+ x^ 

X4 

X3 + X4 
X3+X4 + X5 X3 + X4 



X2+ X3 + X4 
X3 + X4 
X3 + X4 

X2 + X3 + X4 
X3 + X4 

X4 



= X1X5X3 + X1X5X4, 



(65) 



2:3 + X4 



X4 



X4 



X4X5X1 + X4X5X2- 



F = 2{pi.p2)xiX3X5 +2{P2.P3)X2X4X5 + 2{pi.p2 + Pl-Ps + P2-P3)xiX4X5, (66) 

but since 2{pi.p2 + Pi.p3 + P2-P3) = {pi + P2 + P3Y = = 0, 2pi.p2 = {pi + P2Y = s and 2p2-P3 = 
{P2 + Ps)^ = t, we can rewrite F as: 

F = X1X3X5 s + X2X4X5 t. (67) 



3.2.1 Massless case (mi = ... = 777,5 = 0) 

The simplest case is the one that docs not consider mass scales in the diagram. In this situation the 
parametric representation of G has the form: 

00 / X1X3X5 \ ( X2X4X^\ 

The first step is the expansion of the exponentials, after which we find: 

00 

(_\\-D f ni m n2 n2 ni+n2 

G = y „, {sf {tf' / ^ ^ ^ ^ . (69) 

The exponentials are not susceptible to factorization, since the exponents contain only one term. On the 
other hand, in order to apply the method the polynomial U is factorized in the following form: 

U = flf2+flX5+f2X5, (70) 

where we have defined the submultinomials: 
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/i = {xi + X2), 

h = (2^3 + Xi). 

We now proceed to find the multiregion expansion for the factorized multinomial [/, wliich reads: 

1 1 



■S^+ni+n2 



n3,..,n5 

and the expansion of the submultinomials /i and f2'. 



A, 



(71) 



(72) 



(73) 



= (.3 + = ^ ^^^^^^ 4-< V(-n. -n) ' 

118, "9 

where the costraints {A^} are given by the following equations: 

Ai (-J + ni + 11-2 + ^3 + 714 + "5) , 

A2 = (-^3 - ^4 + ng + uy) , (74) 
A3 = (-?^3 - '^s + '^s + "-9) • 
Finally the multiregion expansion for U can be written as follows: 

^ \^ fh ^1 ^2 ^3 ^4 ^5 lij = l 

'^^^a.-.no r(f + ni+n2) r{-m-n,)r{-m-n5)' ^ ' 

then replacing in (j69p and separating the integration variables, we get: 



(76) 



n ■=! r(^,) „,^„, """-'"^ r(f + m + ^2) r(-n3 - n4)r(-n3 - n,) 

J dxi j dx2 y (ix3 y da;4 4^+"^ j dx^ 41+"^+"^+"=, 

and changing the integrations by their equivalents (•), wc finally obtain the prcsolution of G: 

aLir(.,)„^^„/--- r(f +ni+n2)r(-n3-n4)r(-n3-n5)' ^ ^ 

where the constrains {A^} associated to the integrals are: 

A4 = {vi +711+ "e) , 
A5 = {^2 + n2+ 717) , 

Ag = (j/3 + '^i + ns) , (78) 

A7 = (j^4 + '12 + "9) , 

As = (j^5 + 77l + 7l2 + 714 + n^) . 

In this case we generate at most C| = 9 terms or contributions to the solution of G, each of which has 
the form ;-F'(/_i). Nevertheless, three of them do not contribute to the solution, and the remaining terms 
are distributed in two kinematical regions, which wc now summarize. 



17 



Solutions in the region 



< 1: 



G 



Gi + G7 + Gg 



(79) 



Solutions in the region 



< 1: 



G 



— G2 + Gg + Gg 



(80) 



Explicit analytical solution for the kinematical region 



< 1. 



In order to compare our 
s 



< 1. 



results with those of Ref. [TT], we now present the solution which corresponds to the hmit 

The contributions that correspond to the solution of the diagram in this region are associated to the free 
indexes ni, ni and ng. In order to simplify the final resulting expressions we have used the following 
notation Vijk... = Vi + Uj Vk + ■■■■ 

Explicitly the solution in this kinematical region can be written as: 

g(^) =Gi + G7 + G9, (81) 



where the contributions Gi, G7, Gg are given by: 



and the prefactor Xi corresponds to: 



1^1, J^3, 1^12345 
1 + J^1345 - 1 + J/1235 - D 



^ . _Z3 r(i.i2345 - D)T{^ - i/l2)r(f - Z^34)r(f - 1^5)T{D - 1^1345)^(0 - 1, 

Xi ^ [ ^) 



1235 j 



r(i.2)r(i.4)r(i.5)r(^ - i^i2345)r(i? - J^i25)r(i? - 1^345) 



We also have: 



G7 - X7 s^'-'^^'H-'- 3F2 



V2, D- 1^345, -D - t^l45 
1 + D - J^1345, 1 - 1^4 + 1^2 



where the prefactor Xi is: 



_t T N-D r(f - i^l2)r(f - i^34)r(f - V5)T{D - J/I45)r(j^l345 - D)T{V4 - V2) 



V{vi)V{v3)V{v4)V{v^)V{n^ - I/I2345)r(i? - Z^125) 



and finally: 



Gg = Xg s^-^---H-"^ 3F2 I 



D - Vi23,5, 1-V2 + V4 



with: 



X9 = (-1) 



-D 



r(^^i235 - ■P)r(f - ^i2)r(f - ;^34)r(f - v5)t{d - t/235)r(^2 - ^a) 

T{vi)T{u2)T{uz)T(u^)T(^ - i^i2345)r(^ - ^^345) 



(82) 



(83) 



(84) 



(85) 



(86) 



(87) 
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3.2.2 Massive CBox diagram 

The next level of difficulty of this problem corresponds to the addition of another energy scale in the 
topology. In particular let us consider associating masses of magnitude m to some of the propagators of 
the diagram. As a result we now obtain two- variable hypcrgeometric functions as solutions. 

There are different alternatives for factorizing the multinomial U and extending the solution of the 
massless case found previously to a set of massive cases. Nevertheless, not all the possible distributions 
of the mass scale in the propagators generate two- variable series, and in general an arbitrary assignment 
of masses, even if they are equal, produces series solutions of multiplicity fi> 2. The following forms of 
factorization of U allow us to visualize any cases in which it is possible to solve the problem in terms of 
two- variable series. 

Factorization I. 

U = flf2 + flX5 + f2X5, (88) 

where we have defined the submultinomials: 

h = i^^+^^l^ (89) 

With this factorization it is possible to have a two-variable solution for the following massive cases: 
a) 

mi = = m, 

ma = ni4 = = 0. ^ ' 

or 

b) 

mi = m2 m5 = 0, , . 

nis = rrii = m. ^ ' 

Factorization II. Another possible factorization is given by the expression: 

U = f2Xl+ f2X2+ flX5, (92) 

where the submultinomials are given by: 

h = (f^l-^\ (93) 
which allows the evaluation of a more complex case, containing the following mass distribution: 

mi = m2 = 0, .g^. 
ms = 1714 = fnh = m. ^ ' 

Factorization III. We can have yet another factorization of U : 

U = f2X3 + f2X4 + flX5, (95) 

where we have defined: 

./2 = (,/l+X5) 
Jl = [Xl+X2). 

Such a factorization allows to consider the case when we have the following configuration: 
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TO3 = m4 — 0. 



(97) 



Obtaining the presolution in a case with massive propagators. In order to show once again 

the simphcity of NDIM, we will find the multiregion expansion of the diagram G associated to the following 
particular distribution of masses (m\ = m-i = = m) and (ma = 7714 = 0) , and we will show explicitly 
the solution associated to a specific kinematical region. In this case the integral representation in terms 
of Schwinger parameters is given by: 

G = {,-1 ')-■'(-— 

where the multinomial U is given by (|95p : 

U = f2X3+ .f2X4 + flX5, (99) 

with the submultinomials: 

/2 = (/i + X5) and /i = {xi + X2). (100) 
The expansion of the exponentials gives us: 



n-=ir(..)„,^„3 { c/#+-+«3 

where the order in which the multiregion expansions to the multinomials that are present {U, f2 and /i) 
have to be made is simple, and is the following, according to the dependence level between them: 

C/(/2,/l) ^/2(/l) ^/l. 

A little algebra allows to finally obtain the presolution of the diagram G. Explicitly we have that: 

^" n^I>^„,^„,/""--'"" r(f + n, + n,) r(-ni ^ n, ~ n,)n~m - nr) ' ^^^^^ 
where the following constraints have been defined: 



Ai 


= (f + 


n2 + 


"3 + "4 + ?T-5 + "'6 


A2 


= (-rti 


— 714 


- 715 + 717 + "s) , 


A3 


= (-"6 


- nr 


+ 7I9 + 7740) , 


A4 


= Wi + 


112 + 


"9) , 


As 


= {V2 + 


m + 


"10) , 


Ae 


= {V3 + 


112 + 


"4) , 


Ar 


= + 


n3 + 




As 


= (i^5 + 


112 + 


n3 + ne+ n^) . 



(103) 



According to this result we can say that the total number of possible contributions to the solution is 
= 45, and that evidently they correspond to series of multiplicity two. Nevertheless, 18 of these do 
not really contribute due to the particular nature of the linear system build from the constraints. We 
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now show the remaining eontributions, in terms of a set of components Gi^j, according to the type of two- 

p:r:u p:r:u 

variable hypergeometric function (^Fi-^-^ or iT'?-*-" ^ gee appendix) that is generated and also according to 
the arguments of this same function: 



Set. 1 : Gi,5 + Gi.io + Gt.io + Gg.io =^ F^'-^'-^ 



2:1:4 / o Am'^ 

Set. 2 : G5,7 + Gs.g =^ F^'-^'-^ ' ' « ^™ 



s 't 



Set. 3 : Gi,4 + Gi,9 + Gy.g + Gg.g =^ F^ 



1:4:1 / , 4m^ t 



t s 



4,7 + <^4,8 J- ^-^-^ ( |--, — 



Set. 4 : G4 7 + G4 8 ^ 



Set. 5 : G4,5 + G4,io + Gg.g + (Gg.io = 0) ^ f 3;i;i 



1.4.9 / Atti^ t 

Set. 6 : Gi,3 + Gg.r + Gg.S =^ F2:2:0 , — 

\ ' s s 

Set. 7 : Gi,2 + G2,7 + G2,8 =^ F2M \ 

3:2:1 / , a Arn"^ 

Set. 8 : G2,5 + Gs.io =^ F^'-^--^ ' ■ ^ 



4m^ Arm? 



Set. 9 : G3,4 + G30 =^ F4:0:l 



Set. 10 : G2 3 F^ , s t 



It is important to notice that each set of contributions does not necessarily constitutes the final solution 
in certain kinematical region, because it is possible that for some cases the type of solutions get mixed. 
As an example let us consider the case where the magnitude of the variable |t| is the dominant, in which 
case and according to the convergence conditions of each set of solutions, it will be possible to write 
immediately the solution in this region as the algebraic sum of the sets 1, 2, 5, 7 and 8. 

Now let us write explicitly one of the solutions of G. For this purpose we consider the solutions that are 
present in the region where (4m^ > |s|) and (4m^ > \t\), which according to the previous set is composed 
of only one contribution, G2,3: 

g(i&2^1^D=G2,s, (104) 



^ 4to^ ' 4m? 

where the contribution G2 3 refers to the series 



^ (-1)-^ i-ir^+-^ i-mT (sntr ^^'^ 

^'^ UUi'^i'^j) "^Ins! r{D/2 + n2 + n3) r(-ni - 714 - n5)r(-n6 - 7x7) ' ^ ' 

and where the dependent indexes take the following values: 
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ni 
n5 



2^12345 + D 

1/4 



n2 - ris, 



- n3 
= 1^34 — -J, 

"7 = -1^1234 - 

"-9 = -vi - n2 
"10 = V2- ri'i. 



D 



!^345 + f 



- 712 - "3, 
^^2 - "-3, 



(106) 



Making the corresponding replacements and after a little algebra in the expression pOSp . we obtain the 
representation in terms of the Kampe de Feriet generalized hypergeometric function \Fi'- 



KV.U 

':s:v 



G 



t 



Am? ' 4m-^ 
where the prefactor x is given by: 

_Dr(i'l2345 



{a} {a} {c} 

m {-} {-} 



t 



x = {-iy 



D)T{u^23A - f )r(^345 - T)r(T 



Am? ' 4m-^ 



^^34) 



r(i'5)r(i^i25 + 2z^34 - D)V{yi2)T{%] 



(107) 



(108) 



and where the corresponding parameters are: 



Oil 
0L2 



1^12345 


-D, 

D 
2 ' 


/?4 


„ 1 

2 


1^1234 - 


ai 


= >^1, 


^^345 — 


02 


= f^3, 


D 
2 ' 

t'125 1 
2 ' 


.34- f. 


Cl 
C2 


= t^2, 
= J^4, 



2 



J^34 



D 
2 ' 



/3i 

/32 

In this case the series converges if the following condition is satisfied: 

t . 

< 1. 



(109) 



4m2 



Am'^ 



(110) 



3.3 Example III : Four-loop propagator 



In this example the full power of the integration technique here presented is revealed, due to the fact 
that the parametric integral that represents the diagram of Fig. 4 is quite complex. 
The integral representation in momentum space is in this case: 



G 



94 



1 



where the branch momenta h 



ZTT 2 ITT 2 

..,8) are respectively: 



(111) 



Pi, 



ki = qi 
k2 = qi, 
^3 = 92 - 91 
h = 92, 
^5 = 92 - 93 
h = 93- 
fc? = 94 - 93 
h = 94- 



(112) 



The initial parameter matrix associated to the diagram can be easily obtained from its topology: 
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1 



Figure 4: Labelled four-loop propagator. 



M = 



/ Xl 



V 



2:2 + 2:3 
-2:3 




Xl 



X3 



-X3 

■ X4 + X5 
-X5 







2:5 





-2:5 

I- 2^6 4 





Xj 






-Xj 
Xj + Xg 





Xl \ 




Xl j 



The multinomial \J is obtained from the evaluation of the following determinant: 



U = 



U = 



2:1 + 2:2 + 2:3 -3:3 

-X3 X3 + X4 + X5 
-X5 




X5 





-2:5 

-26 + 27 -X7 

— X7 X7 + Xs 



X1X3X5X7 + 2:ia;3a:5a;8 
xixsxexs + X1X4X5XS 



X2X4X5X7 

X2X4X5Xg 

2:22:32:728 ■ 

XiXzX7Xs ■ 

a:22;5a:7a;8 



■ XiX3X7Xs 
X2X4XQX7 
X2X4XQXS 

■ X2X4X7XS 
X3X4X7XS 



■ 2:12:30:62:7 

■ xiX4XeXj 
■ xiX4XeXs 

■ 2:32:40:52:7 
- 2:22:5262:7 

■ 2:22:52:62:3 
2:3x52:6X8 



- X1X4X5X7 - 
■ X2X3X5X8 - 
■ X1X5X6X7 - 
■ X1X4X7X8 - 
■ X3X4X5X8 - 
■ X3X4X6X8 - 
X3X5X7X8, 



■ X2X3X5X7+ 
■ X2X3X6X7+ 
- X2X3X6X8+ 
- X1X5X6X8+ 
- X3X4X6X7+ 
X3X5X6X7+ 



and F is in turn obtained from the evaluation of the determinant of the expression IjllSp : 



F = (X1X2X3X5X7 + X1X2X3X5X8 
X1X2X3X6X8 + X1X2X4X5X8 - 
21X2X3X7X8 + X1X2X4X6X8 - 
X1X3X4X6X7 + X1X2X4X7X8 - 
X1X3X5X6X7 + X1X2X5X7X8 - 
X1X3X5X7X8) pf. 



+ X1X2X3X6X7 + X1X2X4X5X7+ 

h X1X2X4X6X7 + X1X3X4X5X7 + 
h X1X2X5X6X7 + X1X3X4X5X8 + 
I- X1X2X5X6X8 + X1X3X4X6X8 + 
h X1X3X4X7X8 + X1X3X5X6X8 + 



(113) 



(114) 



(115) 



3.3.1 Factorization of the multinomials U and F 

The large number of terms that U{34) and F{21) have makes impossible to get manageable solutions 
if the method is applied without performing the necessary factorizations first. Nevertheless, this type of 
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diagrams that in the massless case are solvable loop by loop, present a form such that they can be easily 
factorized in submultinomials. For this particular case the adequate multimomial factorization is: 



F = xifrpf, 

(116) 

U = xife + fr, 



where the functions fi are given by the equations: 

/? = {X2fe + fb), 

/e = 2:3/4 + (2:4/4 + /s), 

/s = 2:3 (2:4/4 + /a), 

h = 2:5/2 + (.X6/2 + /i), (117) 

/s = 2:5(x6/2 + /l), 

/2 = (xr + xg), 
/i = X7X8. 

3.3.2 Massless case (mi = ... — mg — 0) 

This case is simple to solve loop by loop, but the advantage of the technique NDIM is that it is possible 
to solve the corresponding Feynman integral considering simultaneously all the loops. The parametric 
representation of his diagram is: 

/ Xl.f7 2' 

G = — / — ^ ""'^'^^i ^ (118) 

n'=ir(^.){ {x,h + hV 

where as usual the exponential is expanded, and systematically we make the successive multiregion expan- 
sions associated to the multinomials fi. The order in which these multinomials appear and get expanded 
is: 

u-^h^h^h^h—'h^ h- (119) 

Once the integration process finishes, finally the presolution or multiregion expansion of the diagram G 
is obtained: 

where we have defined the factor: 



r(-ni - n3)r(-n2 - n4)V{-n^ - n7)r(-n6 - n8)r(-n9 - nii)r(-nio - ni2) ' 
and the corresponding constraints are: 



Ai 


-(# + 


ni + 712 + "3) , 


A8 = 


= (i^l + 771 + 7I2) , 


A2 


= (-rii 


- 713 + n4 + «5) , 


A9 = 


= {U2 + 774) , 


A3 


= (-n2 


- 714 + ne + nr) , 


Alo 


= (7/3 + 715 + Tie) , 


A4 


= (-"-5 


- 717 + ns + ng) , 


All 


= {v4 + ns) , 


As 


= {-na 


- 7l8 + Tlio + nil) , 


A12 


= (7^5 + 7lg + 77io) , 


Ae 


= {-ng 


- mi + 7112 + nis) , 


Ai3 


= (z^e + '^12) , 


A7 


= {-nio 


1 - 7112 + 7114 + 7115 ) , 


Ai4 


= {vr + ni3 + 7714) 








Ai5 


= (7^8 + ni3 + 7115) 



(122) 
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The number of possible contributions that the solution has is Cj^ = 1, and it does not correspond to a 
series but to a single term. The solution for this case is simply: 



where the indexes rii are to be replaced by the values: 

m = 2D - Ui - 1/2 - 1^3 - Vi - 1^5 - 1^6 - 1^7 - i^8, = D - U<i - - - I/g, 

n2 = 1/2 + 1^3 + I'A + 1^5 + ^6 + ^^7 + ^8 - 2-D, Uio = Vq + V-j + U^ - D, 

713 = Vi- '^11 = - f , 

Hi = -1/2, ni2 = -1/6, 

715 = ^ ~ V3 ~ 1/4 - 1/5 - l/Q - 1/7 - l^S, 7li3 = Y - V7 - I^S, 

716 = 1/4 + 1/5 + 1/6 + 1/T + 1/s - «14 = 
riT = Z/3 - ^, ni5 =1^7-2-, 
TlS = -1^4- 

3.3.3 A massive case example (mi — ... — niQ — 0) {7717 — ms — 7n) 



(124) 




Figure 5: Mass configuration in the four-loop propagator. 

Let us consider a more complicated situation, which adds one energy scale m to the previous problem, 
as it is shown in Fig. 5. Starting from the parametric representation of the diagram: 

( ^^-2^5 7 exp(-/2(-m2))exp(---^i^p2 
G = ±^ / d-^ ^ ^^^^ + ^^ ^ (125) 

and remembering that /2 = {xr + xs), we proceed to make the expansions, similar to the massless case, 
except that this time an extra expansion appears, due to the exponential that contains the mass scale, 
and which gets expanded later on in the process or replacing the functions At the end of the procedure 
we obtain as presolution the following expression: 

°- fii:;n^..5./" r(D/2+„o «wn,« A„ (129) 
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where is the factor: 



^"■^ r(-ni - n3)T{-n2 - n4,)T{-n5 - nr)T{-n6 - ns)T{-n9 - nii)r(-nio - ni2 - nu) ' 
The new set of constraints {A^}, is given by the following identities: 



Ai = 


(y +ni+n2 +113) , 


Ag = 


= (z/2 + n4> , 


A2 = 


(-ni - ns + 714 + ns) , 


Aio 


= (j^3 + ns + ne) , 


A3 = 


(-n2 - n4 + ne + nr) , 


An 


= {1/4 + ng) , 


A4 = 


(-n5 -nr + ns + ng) , 


A12 


= {1^5 + ng + nio) , 


A5 = 


(-ne - ns + riio + rin) , 


Ai3 


= (j^6 + ni2) , 


A6 = 


(-ng - nil +ni2 +ni3) , 


Ai4 


= (1/7 + ni3 + ni5) 


Ar = 


(-nio - ni2 - ni4 + ni5 + nie) , 


Al5 


= {vs + ni3 + me) 


A8 = 


{ui + ni + n2) . 







(128) 



According to the previous analysis the number of possible contributions to the solution of G is Ci| = 16, 
of which actually only ten contribute to the final solution. These can be distributed in the following 
manner, taking into account the kinematical region of interest: 



Solutions in the region 



Pi 



Aw? 



< 1: 



G 



A. 

Arn? 



Gl + G5 + Gg + G13. 



(129) 



Solutions in the region 



4m 



Pi 



G 



< 1 

Am? 



= G2 + Gq + Gio + Gi4 + G15 + G 



16- 



(130) 



4 Comments 

4.1 The factorization process 

From the point of view of the actual procedure the aspect that optimizes the method has to do with 
how the multinomials that appear in Schwinger's parametric integral are factorized. Is in this process 
where resides the possibility that the integration method be generalized to L loops for both massless and 
massive diagrams. In this respect we worked out a procedure that allows for an adequate factorization of 
the multinomials U and F at the moment in which they are generated, which is something very useful 
since without it the number of sums of these multinomials can be quite large, depending on the number 
of loops and external lines that the topology has. 

In this work we have treated cases in which the multiplicity of the series that conform the solution 
fulfiUes the condition /x = (n — 1), where n is the number of different energy scales present in the diagram. 
With this criterion we have classified in three families the topologies considered here, and for each case 
we give the recipe to factorize the multinomials U and F. 



4.1.1 Diagrams reducible recursively loop by loop 

The topologies that are included in this category are those that in the case in which the theory 
does not contain masses, they can be evaluated loop by loop by successive application of the formula 
corresponding to the one-loop diagram: 
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d^qi 1 



(131) 



whose explicit solution is: 



G 



5(01, a2) 



1 



(132) 



(p2) 



_D : 
2 



where the factor g{ai^ is defined as: 




(133) 



The formulae (|13ip and (|132p can be represented graphically as shown in Fig. 6. 




a-i + 32 -D/2 



Figure 6: Graphic formula for Bubble diagram. 



A situation that happens quite often together with the loop by loop reduction is when there are 
two scalar propagators in series and associated to the same mass scale, with powers ai and 02 in the 
propagators. In this case such powers can be summed and replaced by only one, characterized by the 
power (ai + 02). Graphically this is represented in Fig. 7. 



The massless diagrams which applying formula (|132p can be evaluated recursively loop by loop, are 
also capable of evaluation considering all loops simultaneously and even if different mass scales are added 
to the topology. That this can be done depends only on finding the optimal factorization of U and F. For 
this purpose we use the well known topological analogy between Feynman diagram and resistive electrical 
circuits, in the sense that in both 'something' flows in their branches and 'something' is conserved in their 
vertices. But this analogy goes further, since electrical resistors in an electrical network are equivalent to 
Schwinger's parameters in Feynman diagrams. The analogy is obvious if we observe Fig. 8. 



In this example the following can be seen: in the Feynman diagram one finds that the ratio between 
the multinomials F and U is given by the expression: 



a? + a2 



Figure 7: Reduction of propagators in serie. 




Figure 8: Simple analogy of Feynman diagram and network of resistances. 



27 



F ^ xiX2 ^.^^^ ^ U = xi+X2. (134) 

U X1+X2 

For simplicity we have taken — 1, since the results are independent of this choice. Similarly the 
equivalent resistance of the electrical circuit shown if the figure can be easily evaluated, giving: 

The resemblance of these mathematical structures and also considering the graphic equations shown above, 
allow us to explain the reason why a massless L-loop Feynman diagram can be reduced loop by loop in 
the same way that a resistive circuit can be reduced in terms of sums of series and parallel resistors. Both 
systems have the same topological characteristics. Therefore, using this similitude it is possible to think 
in doing 'sums of parameters' depending on whether we have series or parallel propagators, and with this 
it is possible to find the ratio {F/U) and thus identify directly the multinomials F and U . Nevertheless, 
the main advantage is not really the determination of these multinomials, but the fact that in this way 
it is possible to find simultaneously the adequate factorization for the diagrams solvable directly and 
recursively loop by loop, since in the process of reducing the diagram in terms of series and parallel 
branches, the factorizations of the type {xi + Xj + ... + Xk), adequate for the application of NDIM, are 
generated automatically. This is the trick used in the example (///) in order to find the factorizations. 
A particular case are vacuum fiuctuation diagrams, for which the multinomial F vanishes, and therefore 
for finding U it becomes necessary to assume that the diagram corresponds to a propagator, and for that 
purpose we adequately add two external lines, in such a way that the topology be reducible loop by loop. 
Then the same recipe as above is applied, the ratio {F/U) is determined and finally the already factorized 
multinomial U is found. 

Notice that the factorizations of the type {xi + Xj + ... + Xk) contain the parameters that belong to 
the same loop (once the subtopologics of the loop propagators have been reduced), which seems natural 
since the series-parallel reduction can be done one loop at the time. It is also worth pointing out that 
it is possible to do this 'parameters reduction' independently on whether the propagators are massless 
or massive, which is due to the fact that the mass scales do not appear in the multinomials U and F 
of the parametric integral. Moreover, it is thus clear that when adding equal mass scales, these can 
only be associated to the same loop, and therefore we can be sure that the series solutions increase the 
multiplicity in one unit, which is logical since one more mass scale has been added. As an example, if in 
the multinomials the factor (xi + Xj + ... + Xk) appears, it is possible to add to the topology equal masses 
in various possible situations and associate to the integral of parameters the following exponential factors: 

exp(a;im^), exp((.Ti + Xj) m^), exp((xi + Xj + ... + Xk) m^), etc. 

In the case in which equal masses are distributed in different loops, in most situations it will not be 
possible to factorize adequately the multinomials and the integration method will generate extra sums in 
the series solutions, whose arguments would be one. The way to think about this is that when masses 
are put into different loops this really corresponds to different mass scales and therefore the multiplicity 
of the series solution will grow with the different number of mass scales in the problem. 



4.1.2 Tv^ro loop topologies, with three or more external lines. 

For topologies with more than two external lines, the previous trick is not useful. Nevertheless, there 
is still a generic procedure that can provide an optimal factorization for the multinomial U, which in this 
case is the one that determines and conditions the adequate factorization that must be done in F. This 
procedure can be applied in general to any two-loop diagram, although only in certain topologies optimal 
results are obtained, according to the criteria expressed in ([26|) . In actual practice the factorization of U 
is done directly from the matrix of parameters of the diagram, and since we are considering the two loop 
case, U is obtained from the determinant of the 2x2 symmetric submatrix that relates the loop momenta 
only. The general structure of the determinant is given by the expression: 
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u 



Xa + Xb ±Xb 

±Xb Xb + Xc 



(136) 



where the terms Xa and Xc arc 1-hnear combinations of Schwinger parameters, such that they do not 
contain common parameters among them, that is {^a} H {-^c} = 0- The term Xb is a combination of 
the parameters that are common in the diagonal elements. 

The form of U comes from the determinant of this matrix, that is: 



U = XaXb + XbXc + XaXc, 

where we can do the following factorizations in case they are necessary: 

U={Xa+Xc)Xb + XaXc, 



(137) 



(138) 



U = X,Xb + X,{Xa + X,), 



(139) 



U = Xa {Xb + Xc) + XbXc. (140) 

When adding one more mass scale, the appropriate factorization will depend on the form in which this 
scale gets distributed, and them in the integrand the following factor appears: 

exp(Xam^), exp{{Xa+ Xc)m'^), etc. 

This provides us with a general recipe for the case of two-loop diagrams and more than two external 
points. In our work only three topologies are suitable for applying the method and producing optimal 
solutions, that is they fulfill the condition that the solution multiplicity is less than the number of energy 
scales of the diagram. These are indicated in Fig. 9. 






Figure 9: Two loops diagrams compatible with NDIM. 



4.1.3 The previous topologies + propagator insertions of loop by loop reducible topologies. 

The last case we will consider is a combination of the two-loop topologies shown in Fig. 9, with the 
addition of subgraphs in the loop by loop reducible propagators, and with the advantage that it is also 
possible to include masses in them. In this manner it becomes possible to extend to i-loops the two-loop 
topologies shown previously. 



4.2 Increcise in complexity of the hypergeometric series as the number of 
loops L grows 

In this work we have applied the method to certain L-loop class of diagrams characterized for having 
one, two or three different energy scales, and such that the solutions correspond only to series of zero, one 
or two variables respectively. 
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With this restriction in the type of solution, the topologies to which is possible to apply the technique 
NDIM advantageously are the following: 

• One loop topologies with two, three and four external points, which have been studied in Refs. 

mm- 

• All L-loop topologies of two external points which are loop by loop recursively reducible in the 
massless propagators case, but to which one or two mass scales are added. 

• All L-loop vacuum fluctuation topologies that are reducible loop by loop, and to which up to three 
different energy scales can be added. 

• Some two-loop diagrams, shown in Fig. 9. 



To all these topologies more loops can be added by including one-loop insertions in the propagators, 
and since it is also possible to assign mass to these insertions, we can therefore study a great variety of 
graphs with this integration technique. In conclusion a large family of L-loop diagrams can be evaluated 
with this method. 

Our results indicate that the number of parameters that characterize the generalized hypergeometric 
functions that are obtained as solutions to the diagrams increases linearly with the number of loops L 
of the graph. In the case of topologies that have two energy scales, the solutions are always expressible 
as generalized hypergeometric series of the type ;Li_i, where I corresponds to the hypergeometric order 
and the quantity (21 — 1) is the total number of parameters that it has. Let us see how the number of 
loops L affects the order of the hypergeometric series in the solution. We consider the following examples 
which are associated to the diagram in Fig. 10, and which have two different energy scales {p^,r7i^} and 
L loops. For this class of diagrams we show two topologies which have masses associated to one of the 
loops, with the difference that while the first (A) has one massive propagator the second (B) has both 
propagators massive, and the rest of the propagators are massless. 




Figure 10: Two massive family of L loops propagators. 

The solutions that we have found for these diagrams are given in terms of generalized hypergeometric 
functions iFi-i, where the graphs with the condition (A) have series solutions of the form: 



(2L+1 



(2LH 



2lF(2L-1 
2lF^2L-1) ( ■•• I ^ ) 

while for diagrams of the form (B) we get: 



4m^ 



4mf 
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The example (///) of this work corresponds precisely to a diagram of the type (B), in which L — i 
and whose solutions are just hypergeometric series of type gFg. Another family of solutions that we can 
consider are vacuum fluctuations of the type shown in Fig. 11, composed of two energy scales {m^, M^} 
and which gives rise to solutions of the form: 




Figure 11: Two mass scales of vacuum fluctuations of L loops. 

In case in which there are three energy scales, the solutions can be expressed in terms of double 
hypergeometric series. In the calculations performed in the case of one- loop diagrams 6 , the solutions 
are completely described by Appell functions {Fi, F2, F^, F4), one Horn function H2, and the particular 
cases Si and 5*2 of the de Kampe de Feriet function. Nevertheless, for L > 1 the number of parameters 
that characterize the hypergeometrics grows and then it is not possible to describe the solutions with the 
same functions associated to the 1-loop results and it is necessary to turn to a more general version of the 

p:r:u p:r:u 

two-variable series Fi'-''-^ or F'?-''-^. 

5 Other examples 

Although the previously shown examples are sufficient to demonstrate the versatility and power of 
the method, we want to add briefly other topologies of two, three and four external lines. The last two 
correspond to specific two-loop topologies to which the method is applicable, and which were not discussed 
before. 

5.1 Sunset diagram 



1 




Figure 12: Sunset diagram. 
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U = a;ia;2 + a;ia;3 + a;2a;3. 
F = X1X2X3 pI- 
5.1.1 Massless Ccise (mi = m2 = = 0) 



(-1)- 



r(D/2 + ni)' 



Ai = {D/2 + ni + 712 + na + ^4) 
A2 = {vi + ni + n2 + ns) , 

A3 = {v2 +ni + n2 + rii) , 
A4 = {1^3 + ni + ris + 714) . 



The solution is found immediately, and reads: 



jj T{D/2 - yi)T{D/2 - z/2)r(£'/2 - U3)T{yi ^y^ + y^-p) 

5.1.2 Massive case I (mi = m2 = m) (TO3 = 0) 

^_ (-1)-^ ^^ e^^{-f{-m^))e^v{-F/U) 

For this particular case U has to be factorized as: 

U = X1X2 + 3:3/, 
/ = (a;i + X2) . 



G 



(-1)-^ 



E <l>nu..,n. (— ')"^ (f?)' 



n,ti A, 



r (D/2 + n2) r(-n2 - 714) ' 



Ai = (D/2 + + 713 + rii) , 
A2 = (-711 - 714 + 715 + ne) , 
A3 = {vi + 712 + n3 + ns) , 
A4 = {^2 + 7i2 + n3 + ne) , 

A5 = {^3 +712+ "■4) • 



The solutions, in terms of the contributions Gj, are: 



Gi + G4 + G5 + Gg 
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5.1.3 Massive case II (mi = m2 = to) (tos = M) 



G 



(-1)-^ 



dx 



exp (-/ (-m2) - X3 (-M2)) exp {-F/U) 



U = X1X2 + xaf, 
f = {xi+ X2) ■ 



(149) 



(150) 



G 



(-1)- 



11^=1 r(j^j) 



/ J v^ni,..,n7 



(_^2)"^ (-My (piy^ 



r (£)/2 + na) r(-ni - ng) " 



(151) 



Ai = {D/2 + na + 714 + ng) , 
A2 = (-ni - rig + Tie + 717) , 

A3 = (z/i + ns + n4 + ne) , (152) 

A4 = {1^2 + "-3 + ?T-4 + "-7) , 
A5 — (1/3 + T),2 + W3 + "-5) • 

The solutions determine three kinematical regions, which we write in terms of the contributions Gij: 



' ArrP' \ 2-2-0 

G ( — = Gi,2 + Gi,4 + G2,5 + G2,6 + G2,7 + (G4,5 = 0) + G4,6 + G4,7 =^ FO-sh (...) 



Pi ' Pi 



G 



G 



M 



2 ^2 



Am? ' Aw? ) 



Am pI 



G2,3 + G3,4 =^ F2aa (...) 



(153) 



M2 ' M2 



2:2:0 



, = Gi,3 + G3,5 + G3,6 + G3,7 Folsll (...) 



p:r:u 



In each of these cases the contributions Gij correspond to de Kampe de Feriet functions (F?-*-''). 

5.2 Massless non-plancir Double-Box diagram (On-Shell case) 

U = (a;i + X2) {x5 + xe) + {xa + X4) {x^ + xe) + {xi + X2) {xs + X4) , 

F = xixsxe s + X2X3X5 u + X2X4Xe t, 



(154) 



G = A K,..,r> 



{sr {ur ar 



n"rA, 



T{D/2 + ni + 712 + "3) r(-n4 - n6)r(-n5 - 7i6)r(-7i4 - ng) ' 



(155) 



(-1)-^ 
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Figure 13: Labelled non-planar Box diagram. 



Ai = {D/2 + ni + n2 + + n4 + n5 + no) , 

A2 = (-714 - ne + n7 + ng) , 
A3 = (-^5 - ne + ng + mo) , 

A4 = (-714 - 715 + "11 + ni2) , 
As = {Vl + 771 + 777) , 
Ae = {V2 + 772 + 773 + m) , 
A7 = (j^3 + 771 + 772 + 779) , 
As = (j^4 + 773 + 77lo) , 
Ag = {iy5 + 772 + 77ii) , 

Aio = {ve + ni + ns + ni2) • 

The sets of terms which contribute to the solution are shown as functions of the components Gij and of 
the type of two- valued function that represents them: 



Set 1 : Gi,2 + Gi,i2 + G2,8 + G8,i2 Fi^^ ( 
Set 2 : G2,3 + G2,9 + G3^i2 + G9^i2 Fi-^ 

Set 3 : Gi,3 + Gi,9 + G3,8 + G8,9 =^ F^ 
Set 4 : Gi,ii +G8,ii ^f'^i^o 
Set 5 : G3,ii+G9,ii ^F^ 
Set 6: Gi,io + G8,io =^ F^ ( 



s t 
t ' u 

t s 
s' u 



s u 
u t 



Set 7 : G2 7 =^ 



Set 8 : G7 12 ^ ^^2;o;i 



u t 
t u 



u s 





1 ^ 


t 




' s 


s 


( 


1 ^ 


-) 




' u 


u 
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Set 9 : G2,io = 

Set 10 : Gio,i2 

Set 11 : G3,7 = 

Set 12 : G7,9 = 

Set 13 : G7,io = 

Set 14 : Gio.ii 

Set 15 : G7.11 = 



2:1:2 / 
i;^l:l:0 



1:2:1 / 

^2:0:1 I 



2:1:2 

^1:1:0 



1:2:1 

^2:0:1 



1:2:2 

F2 - - 



2:2 / 
0:0 I 



1:2:2 
F2:0:0 



^1:2:2 



U S 

s u 

" t'7 

t u 
u' s 

u t 

s ' u 
u u 



s s 

t ' u 

t t 

s ' u 



This example has been solved with the usual NDIM method in Ref. TT]. Nevertheless, it is evident 
that doing it the way presented here the solutions are found more directly. They are of the type that 
corresponds to three different energy scales, that is two-variable series, although without generating a big 
number of contributions to the solution, neither redundant sums with unit argument in each one of them. 
Table III allows to visualize more precisely the improvement of our optimization, comparing the solutions 
obtained without any simplification and with the present method. 





Previous NDIM 


New NDIM 


Possible solutions 


6435 


66 


Relevant solutions 


2916 


27 


Multiplicity (/i) 


8 


2 



5.3 Four-loop vertex 




(Table III) 



Figure 14: Labelled four-loop vertex diagram. 
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5.3.1 A case with three energy scales (pf = pi = O) (™6 = Ti? = M) 



G 



dx 



where the functions fi are given by: 



exp (-■ ^^)) exp (-.xg (-Af^)) exp {-F/U) 
F = 2:1X5/2/5 P^, 

u = hhh + hhh + hh, 



/e = XsXg + 0:5/2, 
/5 = X6X7 + hfi-, 

/4 = xe + X7, 

/s = 2:3 + Xi, 

h = xs + xq, 

h = Xl +X2. 



G = 



(-1) 



-2D 



E 



r(f + n3) 



where the factor fl^n} is: 



r(-n4 - n5)r(-n3 -n^- n6)T{-n4 - n6)r(-n3 - 715 - ni2)r(-nio)r(-n4 - nio) ' 



Ai = {D/2 + ns + n4 + ns + ne) , 

A2 = (-714 - 715 + n7 + 778) , 

A3 = (-773 - 775 - 776 + 779 + 77lo) , 

A4 = (-774 - 716 + nn + ni2) , 

A5 = (-7T.3 - n5 - 7712 + 7713 + ^4) 
A6 = (-7710 + ni5 + 7716) , 
A7 = (-774 - Tlio + 7li7 + nis) , 
As = {1^1 +713+ "7) , 



Ag = (7^2 + 778) , 

Aio = (i^3 + 7715) , 

All = (j^4 + 7716) , 

A12 = Ws + n^ + nu) , 
Ai3 = {1^6 + 771 + ng + 7717) , 

Ai4 = (j^7 + 779 + 7718) , 

Ai5 = + nil + nia) , 
A16 = (J^9 + n2+ nil + ni4) . 



The sets of solutions are given by the foUowing equations in terms of the contributions Gij. 



Set 1 : Gl,2 + (jl,5 + Gl,13 + G2,4 + G2,18 + (G4,5 = 0) + G4,13 + G5,18 + Gl3, 
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3:3:2 



3:3:2 



2 2 
771 p 



Set 2 : Gi,3 + G3,4 + ^3,18 =^ F2.Z.2 { | ^, 
Set 3 : G2,3 + 63,5 + 63,13 =^ F^^M ' 1 ^ 



2:4:3 



Set 4 : 6*7,9 + ^9,12 =^ F3\2\i 



( 1 


M2\ 







Set 5 : G7,ii + (Gii,i2 = 0) ^ F4 



,3:3:3 



M2' p2 
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^2:4:3 



Set 7 : + G4,ii =^ F^''^'-^^ 



Set 8 : Gil is =^ F^'-^^'^ 



( 1— , 






M2 j 



/ 1 






M2 J 



f 1 — 


M2\ 




p2 j 



—3:2:4 / , 

Set 9 : G2,9 + G5,9 ^ F2:2:2 _ 
Set 10 : G9,i3 =^ F^"^''-^ 



Set 11 : Gi,7 + Gi,i2 + €4,7 + G4,i2 =^ i^^l;? 



_2:3:3 / , AfS 2 \ 

Set 12 : G7,i8 + Gi2,i8 =^ F^-^-^ I \ - -!!^\ 



( 1 "^^ 


M2\ 







p2 ' ^/2 y 



Set 13 : G3,9 ^ F2^2:2 -^^--T 



Set 14 : G2,7 + G2,i2 + Gg.r + G5,i2 ^ i^^;!;? 



Set 15 : G7,i3 + Gi2,i3 =^ F^^^ 





1 M2 


m2\ 


( 


9" ' 

















_4:2:3 / -2 ^2 

Set 16 : G3 11 ^ F3^2:i 



I p m 

I ^'~M2 ; ■ 

In this example the fact that (M^ > m^) or vicevcrsa leaves out automatically sets of solutions, since 
only one of these conditions can appear in the same topology. 



6 Conclusions 

The integration method here presented is a powerful tool for solving Feynman integrals with arbitrary 
exponents in the propagators. Its main characteristics resides in its simplicity in finding the solution of 
diagrams, even in the evaluation of L-loop diagrams. This technique, which can be more appropriately 
called Integration by Fractional Expansion, is clearly an important competitor to other Feynman diagram 
evaluation methods. In spite of its generality, we have found that the technique is adequately applicable to 
certain types of graphs, obtaining solutions which we call optimal, in the sense that they are algebraically 
manageable, and whose complexity does not depend on the number of loops L of the diagram. Although 
this method can be applied to any diagram, the complexity of the solutions in the more general cases 
has a high dependence on L, which implies a large number of terms as well as a high multiplicity of 
the resulting series expansions, which makes the analysis of the solution highly non trivial. Imposing to 
evaluate only topologies whose solutions can be represented as series in which the variables are expressed 
as ratios between the different energy scales of the diagram, determines that the multiplicity of each of 
them fulfills the relation jx = (n — 1), where n is the number of different energy scales in the diagram. 
This is the restriction that allowed us to find that the method here proposed is applicable in a natural 
way to certain L-loop class of diagrams, which in the massless case happen to be solvable loop by loop. 
Using this condition it becomes possible, depending on the number of energy scales present in this class of 
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topologies, obtain solutions as expansions in multivalued series. Nevertheless, we have considered adequate 
to evaluate topologies whose solutions could be expressed at most in terms of two-variable series, because 
hypergeometric series of this type have boon amply treated in the literature. The advantages of this 
integration method can be summarized in the following: 

• For the family of L loop graphs, reducible recursively loop by loop, the application of the technique 
makes possible to find the solution considering all loops simultaneously, through its Schwinger 
parametric representation. 

• The integration method allows to extend in a simple manner the solutions to cases in which different 
mass scales are considered in the graph. 

• It is possible to apply the integration method to certain three and four external line topologies, and 
generalize it to L loops making one- loop insertions in the propagators. In general the application of 
this technique can be extended to more external lines, but the solutions are series with > 2, and 
therefore have not been considered in this work. 

• Independently of the number of loops, in topologies of the loop by loop reducible type, the mul- 
tiplicity /i of the hypergeometric scries solution docs not have any dependence on the number of 
loops L, and only depends on the number of different energy scales that characterize the graph. For 
topologies that are not in the class reducible loop by loop, except for the special cases described 
previously, the multiplicity grows when L increases. This fact makes this type of topologies not op- 
timal for the application of the method, because the complexity of the solutions grows very quickly 
with L. We can estimate this growth in terms of the contributions that the solution contains, using 
the following expression: 



Number of contributions of the solution ~ ■^rrr^^TTTTr" • 

A significant fact is that the method allows to easily add masses to any type of topologies. In general 
for generic diagrams we have been able to conclude that: 

• If we add M different mass scales, the multiplicity of the series expressions that are part of the 
solution grows in M units. 

• In most cases if we add masses to propagators belonging to I different loops, the masses have to be 
considered different, even if they have the same value, and the multiplicity increases in I units. 

• The conjunction of the previous premises implies that if one associates M equal masses to the 
propagators of a loop, the multiplicity of the series solution only increases in one unit. 
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A Mathematical formalism of NDIM 



In section (2.2) the technique was introduced, and also its origins starting from the identity associated 
to the loop integral parametrization and known as Schwinger's paranietrization: 

oo 

-^^Y^Jdx x^-i eM-Ax), (161) 



and furthermore the equivalence between the integral symbol and a delta Kronecker psp was shown 
starting from the equation (|16ip : 

dx x^+--' EE r (/3) ^pYy^ Sp+n,o. (162) 

For simplicity we have eliminated the limits of the integral, since this identity has validity only in the 
context of the expansion of the integrand in p6ip . This expression is crucial for the development of the 
method NDIM, since in the evaluation of Feynman diagrams the corresponding Schwinger parametric 
representation is a generalized structure of the expression (I161|) . 

A.l Some properties 

In order to study some of the properties of (|162p . it is convenient to introduce a useful notation for 
helping us to formalize the mechanism of the technique NDIM. Thus: 

dx x'^i+'^^-i = (i/i + 1^2) , (163) 
where vi and 1^2 are indexes that can assume arbitrary values. 

A. 1.1 Property 1. Commutativity of indexes 

We can write (I163P explicitly in two possible forms, according to formula (|162p . and then: 

{vi+V2)= { (164) 
, r(^i + l) ^ 

Starting from (|16ip it is possible to show the equivalence of both forms in ()164p , and for this purpose it 
is enough to expand the exponential of the integrand and replace (•) for this case. 
Let us consider the following integral representation: 

00 

^ " / "^"^ ^""P ^"^"^^ ' ^^^^^ 



which in terms of series, in the sense proposed (|163p . turns out to be: 

3:7 = FMi:Ffe^-^" + (166) 

F {^y ) F (y 1) 

Selecting now {vi + V2) = — ^it^ <^ 1^1+1^2,01 directly obtain the equality in ([165]). Analogously 

F (v ) F (v -j- 1) 

selecting {vi + 1^2) = — -^yi 1^1+1/2,0 and using afterwards the identity: 
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r(.) ^,_j,-.r|i+^_ 1^^^^ 



r{y-z) r(i-s) 

with y = vi and z = 2vi, the equahty (|165p is obtained once again, and therefore the equivalence between 
both forms of writing {vi + V2) is quickly established. 

Another form of writing p63p . which is quite useful for simplifying terms that contain the factor: 

(168) 



r(m+ 1)' 

where m is an arbitrary index, is the following: 

Tim + 1) 

{vi + V2) ^ {vi + V2 - ra + m) = (-to + to) = T{-m)———— S^^+^^^q. (169) 

(-1)™ 

Notice that we have made explicit use of the Kronecker delta in order to simplify the parenthesis (•), and 
since this is in the context of expansions, the Kronecker delta remains as an indication of the constraints 
between the indexes vi and V2. 



A. 1.2 Property 11. Multiregion Expansion (MRE) 

Consider the binomial expression: 

{A,+A2f\ (170) 

where the quantities Ai,A2 and can take arbitrary values. In this case there are two limits or possible 
regions with respect to the expansion: the region where (^1 > A2) and the region where {Ai < ^2)- These 
expansions are respectively: 

1. Region {Ai > A2) 



iA, + A2f'' ^Af'^Y. 



71=0 ^ 



1) 



A2 
A, 



(171) 



2. Region {Ai < A2) 



n— ^ ^ ^ 

The factor {u)n is called Pochhammer symbol and is defined as: 

= (173) 

We have thus obtained separate expansions in the two possible limits. Nevertheless, it is possible to 
express both results using a single series that contains both limiting regions. In this sense we can say that 
this type of expansion corresponds to a multiregion series representation of the binomial (|170p . through 
the use of the integral representation of the denominator indicated in (jl6ip . Then we have: 

00 

{Ai + Aa)^" = -^^^ J dx x'f"-^ exp(-x^i) exp(-a;A2), (174) 


and the exponentials are expanded separately, with the result: 
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(_^ \ni+n2 

r(ni + l)r(n2 + 1) 1 ' 



Making use of the identity ()163p we get the muhiregion binomial expansion: 



(175) 



(176) 



where the parenthesis (•), according to property (|164p . can be expressed this time in three different ways, 
ahhough in each one of them we wiU have the same Kronecker deha which ehminates one of the two sums 
that are present: 



(=Fi^ + rii + 712) = < 



r {n2 + ni) 



r(n2 


+ 1) 


(-1 


^712 


r(ni 


+ 1) 


(-1 






+ 1) 


(-1) 





(177) 



delta, can be found evaluating the combinatorics C^gJ^as' which in this case is = 2. Let us see what 
happens if one sums with respect to a particular index: 

1. Sum respect to n2 

Using for this case the following equality: 



(-1)' 



and then replacing in ()176|) . we get: 



(?ll + 1) 



or equivalently: 



ni— ^ 



(178) 



(179) 



(180) 



which gives the expression associated to the region {Ai < A2), obtained previously in (|172p . 
2. Sum respect to ni 
Analogously, we now use the identity: 



and replacing in p76p . one gets: 



- r(7i2 + 1) V ^1 



A2^"= 



(181) 



(182) 
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an expression that was already found in (|17ip . and valid in the region {Ai > A2). 

The fundamental idea that has been exposed in the previous demonstration is that using the definition 
(jl62p one can do a binomial expansion that differs from the conventional expansion, in the sense that one 
obtains an expression where all the possible limits of the binomial are included simultaneously. 

A. 2 Multiregion expansion of a multinomial 

The previous discussion can be generalized to the multiregion expansion of a multinomial of ni terms: 

(A, + ... + ^ E - E ^nu..^. Ar...Ar ^^''^;i^-^''^\ im 

where for simplicity we have defined the following notation: 

(h = (_l)ni + ...+ni 1 Qg^N 

^"i--"' ^ ' r(ni + l)..r(n, + l) ^ ' 

The number of expansions that can be obtained starting from (|183p is given at most by all the possible 
forms of evaluating some of the sums, using for this the Kronecker delta generated by the same expansion, 
that is C"' = n; possible forms. More generally, a function expressed as a multiregion expansion through 
a sums and S Kronecker deltas, can be evaluated at most in: 

possible forms, each one of them expressed in terms of series of multiplicity /z = {a — 5). All the resulting 
series are representations with respect to the ratios between the terms of the multinomial, and all of them 
correspond to multivariable generalizations of the hypergeometric function. 



B Negative dimension D ? 

The original name of the integration method here presented, NDIM (Negative Dimension Integration 
Method), comes from applying the expansion and subsequent association, integral<^Kronecker delta, over 
the gaussian integral in D dimensions: 

exp(afc2), (186) 

in 2 

where k corresponds to a 4-momentum and a is an arbitrary parameter. This is a typical integral that 
appears upon using Schwinger's parametrization and after using the completion of squares procedure of 
the loop momenta. Starting from this equation it is possible to deduce in a similar way to equation (jlSp 
the following identity ^ : 

^(fc2)"=n!<5„+.,o. (187) 

Since this is in fact a Taylor expansion, we assume that n > 0, and therefore the constraint associated to 
the Kronecker delta in (|187p requires that < 0, that is, la dimension has to be negative, which is then 
the origin of the name of this integration technique. Nevertheless, there can appear reasonable doubts 
with respect to some concepts that the NDIM uses in its original approach. 

It is entirely possible for D to be a negative integer considering that the Feynman integrals are analytic 
functions in D arbitrary dimensions, but rigorously what is wanted is to solve the loop integrals in the 
limit D +4 and not in the limit D — > —4. The inconsistency is even more clear if the dimensional 
regularization prescription is used, and the dimension D now includes a non integer piece, the dimensional 
regulator e, which cannot be associated to the sum index n in (jl87p . 
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Such inconsistencies could be resolved if the expansion of the exponential could be written as: 



, , , x"^ , , d"exp(a;) 

expfx) — lim > — such that 

^ ' C^o r(m + l) dx™ 



= 1, (188) 

x=0 



where |^| < 1 and m is an index that increases in one unit steps. Such an expansion is possible, but 
for this purpose it is necessary to resort to the area of calculus called Fractional Calculus [531 Ull US] , 
through which it can be shown that both the derivative and the integration operators can be represented 
in terms of a unique operator, and in second term shows that the analytical continuation in the operator 
order can in fact be done, and can be fractional and even complex. This is very useful in order to justify 
the validity of the NDIM integration technique, and therefore it is necessary to review some concepts for 
understanding this type of calculus. 



B.l Preliminaries : The fractional expansion of the exponential function. 



We are interested in knowing the nature of the expansion that is performed in the exponential 
function in (|186p . since due to the arguments mentioned before it should not really correspond to a Taylor 
expansion. 

In order to provide a more rigorous basis for equation ()188p . it becomes necessary to define the 
integration and derivative operators in a generalized form. In fractional calculus is possible to define the 
fractional integral of order a as follows: 



^D-Vix) 



1 



Via) 



dt 



fit) 



{x-ty- 



(189) 



where c G M and a is an arbitrary quantity. For the particular case in which c = in ()189p . we get the so 
called Riemann-Liouville fractional integral: 



oi5-«/(x) 



T{a) 



dt 



fit) 



ix-t)^ 



(190) 



Another version of the fractional integral can be obtained making c = ~oo, which is called Liouville 
fractional integral: 



^D-Vix) = 



1 



r(a) 



dt 



fit) 



(x-t)i 



(191) 



This last version of the fractional integral will be used in this work in order to justify more rigorously 
the integration method of Feynman diagrams here presented. Certainly this would be incomplete if we 
do not define also the fractional derivative (in the Liouville version). For an order a which fulfills that 
(m — 1) < a < m, with m £ N, the fractional derivative is defined by the expression ^aoD" f{x) — 

[_oo-D™~"/(x)], or in operational terms: 



1 



.D:fix) = 



r(m — a) dx" 



dt 



fit) 



(x-t)^ 



(m — 1) < a < m, 



(192) 



dx^ 



:fi^) 



We are particularly interested in the Liouville integro-differcntial operator, so for simplifying the notation 
let us define _ooZ?" = D", where we can have (a < 0), which is an integral, or (a > 0), a derivative. 
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A property that makes this version of the fractional operator particularly interesting is the effect that 
it has over the exponential function, and which in practice is an extension of what is done by integrals 
and derivatives of integer order (conventional calculus) over this function: 

exp{(3x) a < (integration of order a), 

(193) 

/?"exp(/3x) a > (derivation of order a). 
This expression can be easily proved using the equations ()191|) and p92|) . 

B.2 Taylor and Taylor-Riemann series. The fractional expansion of diagram 
G 

Particular interest has the series representation of the exponential function, since it is the only 
function over which the expansions in the Schwinger parameter integral are done, either because it appears 
explicitly in the integrand, or because it is used implicitly in order to do the multiregion expansions of the 
multinomials present in the problem. When the expansions in the integral of the parametric integral are 
done, we automatically assume that this is a Taylor series. Nevertheless, what is really done implicitly is 
a Taylor-Riemann expansion of the exponential, which reads: 

exp(a;)= V [exp(a;)] I , — ^- -. (194) 



n — — oo 



The parameter ^ is arbitrary and fulfills the condition < |^| < 1. Notice that taking ^ ^ the Taylor 
series for the exponential is obtained. Moreover, since D"^^ [exp(a;)] = 1, according to formula (|193p . 
one gets: 

exp(a;) = > — -. (195) 



n— — OQ 



If we make now the change of variables to = ti + ^, we can rewrite the previous expansion as: 

exp(.)= E (196) 

which looks like the conventional Taylor series expansion of the exponential, with the difference that the 
index of the sum has been analytically continued to negative and fractional values. This justifies that both 
D and the propagator powers {i^i, ...,^Ar} do not change its nature due to some analytical continuation 
required by the Kronecker deltas generated in the process, but rather are the expansion indexes the ones 
that can acquire arbitrary values, including negative and fractional. But moreover it justifies the presence 
of Gamma functions in the multiregion expansion of the diagram, both in the numerator and denominator, 
which having the form V[—n), where n is an index associated to a sum, give a multiregion expansion not 
necessarily infinite or vanishing. 

After having said that, it is still necessary to justify what refers to the integration method, since they 
correspond to hypergeometric series whose summation indexes are positive integers. This can be explained 
by saying that for applying the integration method NDIM an analytical continuation of the summation 
indexes has to be done, or equivalently the exponentials are expanded in Taylor-Riemann series, using the 
expression (|196|) . In practical terms it is enough to assume a fractional arbitrary parameter ^ associated 
to the indexes of all the sums, and at the end of the integration process, when each contribution to the 
solution is determined, evaluate the limit: 

oo-l-^ oo 

rn— — oo+^ m— 
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which naturally happens since all the solutions that are obtained have the factor — -, and 

Tirti +1 + 4) 

the limit ^ — > we have that: 



^ I — , para m > 0, 

lim \ , = { (198) 
Q , para m < 0. 



5^0 r(m + 1 + C) 



Therefore, strictly speaking we should write the Multiregion Expansion (|2ip of a generic diagram G in 
the following form: 



G= (-i)-^iim5^,..,^^^o E - E <^n....,n„ n(Q,')"\ n (-m,^^- 



7 — 1 1 — -I- 1 



(199) 



(z^,-+«,-} 1^ (/3,-+7,) 

i\ r(z.,) i\ r(/3^.) ■ 



Operationally NDIM does not consider the presence of the fractional part ^ in the summation index nor 
its extension to negative values. Nevertheless, the analytic continuation of the summation indexes and 
the evaluation of the limit at the end of the process, are inherent to the process and justify why it works. 
A confusion might arise from the similarity of the Taylor series expansion of the exponential and its 
Taylor-Riemann expansion (equation (jl96p ). something that explains why both versions of the expansion 
give the same result, which in turn shows that there is a conceptual but not operational difference. Thus 
it would be more appropriate to call the method integration by fractional expansion (IBFE) instead of 
NDIM. 



C Hypergeometric Functions of one and two variables 

In this work we have developed the evaluation of loop integrals of certain class of diagrams with 
arbitrary L. All the found solutions have been presented in terms of hypergeometric series of one or 
two variables, which correspond in a natural way to an serie representation with respect to the ratio 
between the energy scales associated to the graph. The purpose of this appendix is precisely to provide 
the necessary information about the hypergeometric functions [26l [27l [28l [29] and especially with respect 
to their convergence properties. 



C.l Definition of the generalized hypergeometric function 

In those cases in which the solutions contain two energy scales the series representations are expressed 
in terms of generalized hypergeometric functions: 



(ai, a^; 6i, 6g_i; z) = g-Fg-i 



{«} 
{6} 



k 



V (ai)fc-(a9)fc ^ ,200) 
^^^(M,...(6,-r),fc!' ^'"'^ 



where the factors (a)^ are called Pochhammer symbols, defined as: 

The convergence conditions for these series are given according to the magnitude of the argument: 

1. If |z| < 1 the series converges absolutely. The variable z represents the ratio between energy scales 
of the topology. 



45 



2. If = 1, the necessary requirement for the convergence of the series is that Re (w) > 0, where w is 
called parametric excess and is given by: 



q 9+1 

For the convergence when = — 1 it is sufficient that Re (w) > —1. 



(202) 



C.2 Some identities of the Pochhammer symbols 

The following identities are useful when building the hypergeometric function from the contributions 
that in turn are obtained from the multiregion expansion of an arbitrary diagram G. The Pochhammer 
symbols are defined as follows: 



n(«+i) 

3=0 
1 



si n > 0, 



(203) 



SI n 



0, 



Nevertheless, the series that are the result of applying the NDIM to Feynman integrals will always contain 
factors of the type r(o ± n) and T{a ± 2n). In these cases it is convenient to use the following formulae 
for finding the hypergeometric representations: 



(«)-n = 



r(a- 



r(a + n) 

-n) {-ly 



T{a) 



(1 - a)n ' 



(a) 



2n 



r(a + 2n) 

no) 



(204) 
(205) 
(206) 



C.3 Two- variable hypergeometric functions 

In this part of the appendix we describe the functions which are useful for representing those cases 
that consider three different energy scales, in the general case L > 1. Such functions correspond to the de 

p:r:u 

Kampe de Feriet generalized double hypergeometric function Fi-^-'" and the generalized hypergeometric 

p:r:u 

pq:s:v_ rj:^^ gj.^^ defined as: 



pqfs^l ai,..,a,. ci,..,c„ 

bi,..,bs di,..,dy 



x,y 



' m {b} {d} 



x,y 



u 



n i<^j)n+m n {aj)n U fe)™ „ ^ 



(207) 



3=1 3=1 j=l 



n! ml 



where the convergence of the double series requires that the following relation between the indexes is 
fulfilled: 



p + r^q + s + 1, 



(208) 
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p + u^q + v + 1, 
and also that the arguments fulfill the condition: 

1 1 

|a;| + \y\ < 1 , for the case in which {p > q), 

max{\x\ , |y|} < 1 , in the case {p ^ q). 

For the other series that occurs frequently we have the following definition: 



(209) 



'^•.s:v ( O^ij-.^Ckp dij-.-jdr Ci,..,C^i 



PlT-^Pq bi,..,bs di,..,dv 



x,y 



_^-^( W {a} {c} 



m {b} {d} 



x,y 



n—m 



i=i 



a;" y"* 



n! m! 



(210) 



where the convergence of the double series requires that the following relations between the indexes are 
satisfied: 



p + r^q + s + 1, 



(211) 



q + u^p + v + 1. 



(212) 



All solutions obtained with the technique of integration NDIM fuUfiU the previous condition in the sum 
indexes. The conditions of convergence of arguments may be worked out using Horns general theory of 
convergence [30] , 



47 



D Bibliografia 
References 

[1] V.A.Smirnov, Evaluating Feynman Integrals (Springer, Berlin, Heidelberg, 2004); and references 
therein. 

S.Moch, P.Uwer, S.Weinzierl, J.Math.Phys. 43 (2002) 3363, (|hep-ph/0110083). 



T.Huber, D.Maitre, Comput.Phys.Commun. 175 (2006) 122-144, ( |hep-pli/0507094) . 
M.Yu.Kalmykov JHEP 0604 (2006) 056, (| hep-tli/0602"028| . 
I.G.Halliday, R.M. Ricotta, Phys. Lett. B 193 (1987) 241. 



C.Anastasiou, E.W.N.Glover, C.Oleari, Nucl. Phys. B 572 (2000) 307, ( ,hep-ph/9907494[ ). 



A.T.Suzuki, E.S.Santos, A.G.M.Schmidt, Eur.Phys.J. C 26 (2002) 125, (hep-th/0205158|. 



A.T.Suzuki, E.S.Santos, A.G.M.Schmidt, J.Phys. A36 (2003) 4465, ( hep-ph/0210148[ ). 
A.T.Suzuki, A.G.M.Schmidt, J.Phys. A31 (1998) 8023. 



C.Anastasiou, E.W.N.Glover, C.Oleari, Nucl.Phys. B 565 (2000) 445, ( |hep-ph/9907523| . 

C.Anastasiou, E.W.N.Glover, C.01eari,Nucl.Phys. B 575 (2000) 416. [Erratum-ibid. B 585 (2000) 
763]. 



A.T.Suzuki, A.G.M.Schmidt, Can.J.Phys. 78 (2000) 769, (hep-th/9904195|. 



A.T.Suzuki, A.G.M.Schmidt, JHEP 9709 (1997) 002, ( |hep-th/9709024D . 



A.T.Suzuki, A.G.M.Schmidt, Eur.Phys.J. C 5 (1998) 175, ( |hep-th/9709144| . 

A.T.Suzuki, A.G.M.Schmidt, Solutions for a massless off-shell two-loop three-point vertex, 
( |hep-th/97T2104| . 

A.T.Suzuki, A.G.M. Schmidt, Phys.Rev. D58 (1998) 047701, (|hep-th/9712"]M). 



A.T.Suzuki, A.G.M.Schmidt, J.Phys. A35 (2002) 151, (hep-th/0110047). 



I.Gonzalez, I.Schmidt, Phys.Rev. D72 (2005) 106006, ( |hep-th/0508013| . 
E.E.Boos, A.I.Davydychev, Theor. Math. Phys. 89 (1991) 1052. 
A.I.Davydychev, J. Math.Phys.32 (1991) 1052. 
A.I.Davydychev, J. Math.Phys.33 (1992) 358. 



A.T.Suzuki, A.G.M. Schmidt, J.Comput.Phys. 168 (2001) 207, ( hep-th/0008143D . 

R.Hilfer, Applications of Fractional Calculus in Physics (World Scientific Publishing Co. Pte. Ltd., 
Singapore, 2000). 

K.S.Miller, B.Ross, An introduction to the fractional calculus and fractional differential equations 
(John Wiley & Sons, Inc., 1993). 

K.B.Oldham, J.Spanier, The Fractional Calculus, Theory and Applications of Differentiation and 
Integration to Arbitrary Order (Academic Press, Inc., 1974). 

W.N.Bailey, Generalized Hypergeometric Functions (Cambridge University Press, Cambridge, 1966). 



48 



[27] G. Gasper, M.Rahman, Basic Hypergeometric Series (Cambridge University Press, 1990). 

[28] L.J. Slater, Generalized Hypergeometric Functions (Cambridge University Press, 1966). 

[29] L.S.Gradshtcyn, L.M.Ryzhik, Table of Integrals, Series, and Products, Sixth Edition (Academic 
Press, 2000). 

[30] H. Exton, Multiple Hypergeometric Functions and Applications (Ellis Horwood,Westergate, England, 
1976). 



49 



